Measurement of the decay of Fock states in a superconducting quantum circuit 
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We demonstrate the controlled generation of Fock states with up to 15 photons in a microwave 
coplanar waveguide resonator coupled to a superconducting phase qubit. The subsequent decay of 
the Fock states, due to dissipation, is then monitored by varying the time delay between preparing 
the state and performing a number-state analysis. We find that the decay dynamics can be described 
by a master equation where the lifetime of the n- photon Fock state scales as 1/n, in agreement with 
theory. We have also generated a coherent state in the microwave resonator, and monitored its decay 
process. We demonstrate that the coherent state maintains a Poisson distribution as it decays, with 
an average photon number that decreases with the same characteristic decay time as the one-photon 
Fock state. 
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Fock states are the eigenstates of the harmonic oscil- 
lator, and as such play a fundamental role in quantum 
communication and information. Due to the linear na- 
ture of the harmonic oscillator, transitions between adja- 
cent Fock states are degenerate, making it impossible to 
generate and detect such states with classical linear tech- 
niques. Interposing a non-linear quantum system, such 
as that provided by a two- level qubit, allows such states 
to be probed unambiguously 0, H & |, S @, 0, ! S M 
Hi 12, ID, 11, 15 1 . However, the deterministic generation 
and characterization of Fock states with arbitrary pho- 
ton numbers has not been demonstrated until recently 
Q, using a superconducting phase qubit coupled to a 
microwave resonator. The exquisite control provided by 
this implementation allows new tests of the fundamental 
physics of Fock states. 

In this Letter, we present the first experimental mea- 
surement of the time decay of Fock states in a super- 
conducting circuit. Although measurements of the relax- 
ation dynamics have been done for the 2-phonon Fock 
state of a trapped ion coupled to a zero-temperature am- 
plitude reservoir [l6j |. and proposed for higher n Fock 
states of photons by performing a statistical analysis of a 
sequence of quantum non-demolition measurements Q, 
a direct measurement beyond n — 2 has not been at- 
tempted to date. The experiment described here is per- 
formed with a system similar to that in Ref . 0] , but using 
a microwave resonator with significantly better perfor- 
mance, allowing the generation of Fock states with up to 
15 photons. The subsequent decay of the Fock states, as 
well as the decay of a coherent state generated by a clas- 
sical microwave pulse, have been directly measured and 
shown to display the expected dynamics. More impor- 
tantly, we confirm that the lifetime of an n-photon Fock 
state accurately scales as T n = T\/n, where T\ is the life- 
time for the n = 1 Fock state [17| . The lifetime T n scales 
as 1 /n because each of the n photons can independently 
decay. 



In our experiment, the Fock states are generated by 
first driving the phase qubit from the ground \g) to the 
excited |e) state, during which the qubit and resonator 
are tuned out of resonance. We then turn on the coupling 
between the qubit and resonator by frequency-tuning 
them into resonance, allowing the qubit to transfer the 
photon to the resonator, generating an n = 1 photon 
Fock state in the resonator and leaving the qubit in its \g) 
state. The qubit is then de-tuned from the resonator, so 
that the subsequent photon dynamics are determined by 
the resonator parameters alone. Because this sequence 
can be repeated n times to produce an n-photon Fock 
state, this method is scalable to arbitrary photon num- 
bers, limited only by the resonator decoherence time, the 
photon transfer speed through the qubit-resonator cou- 
pling, and the fidelities of the qubit |e) state preparation 
and the qubit-resonator tuning pulse. A resonator with 
a relatively long T\ is thus desirable for both preparing 
high n Fock states and measuring the subsequent decay 
dynamics. 

Our device is composed of a half-wavelength coplanar 
waveguide resonator, capacitively coupled to a supercon- 
ducting phase qubit, as shown in Fig. [TJ It has a layout 
similar to that used in previous experiments 0, EH • In 
order to improve the resonator lifetime T±, we have im- 
plemented a new design and fabrication procedure, that 
reduces the amount of amorphous dielectric in the vicin- 
ity of the resonator, thereby reducing microwave dielec- 
tric loss [l9j]- We also used a resonator fabricated from 
superconducting Re (with T c « 1.8 K), as this metal 
has substantially less native oxide than Al, which was 
used previously. As depicted in Fig. [U the fabrication 
incorporates two resonator designs on the same wafer, 
allowing us to measure the resonator quality factor Q 
from a classical resonance measurement [l^|, as well as 
the resonator T\ using a qubit-based measurement 0]. 
The layout of the two resonator designs are as identical 
as possible, to match any dissipation due to (unwanted) 
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FIG. 1: Photomicrograph of a coplanar waveguide resonator 
coupled to a phase qubit. The total length of the half- 
wavelength resonator is 8.8 mm. The insets show the two 
designs of the microwave lines that were varied in fabrica- 
tion to connect the resonator either to a qubit or a classical 
preamplifier, allowing qubit-based Ti or classical resonator Q 
measurements, respectively. 



radiating modes in the microwave circuit. 

The qubit is characterized using pulse sequences de- 
scribed previously [2(j. We found an energy relaxation 
time 600 ns, a phase coherence time « 100 ns, and 
a measurement visibility between the qubit \g) and |e) 
states of close to 90%. The initial characterization of the 
microwave resonator was made using a pulse sequence 
that swaps a photon into and out of the resonator, as de- 
scribed above and in Ref. 



2l| . The decay of the resonator 



state measured in this manner is plotted in Fig. EIJa), 
showing a decay time T[ = 3.47 /is that is three times 
longer than in our previous experiment Q- The decay 
time for a Ramsey sequence gives T' 2 = 6.53 fj,s, approx- 
imately twice the resonator T{, which indicates a long 
intrinsic dephasing time for the resonator, 3> Tj, as 
expected. In addition, the classical measurement [19( of 
the resonator quality factor Q ss 130, 000 at low power 
predicts T\ = Q/oj r — 3.15 /xs, consistent with the qubit 
Ti measurement. Here, u r /2ii — 6.570 GHz is the res- 
onator oscillation frequency, in agreement with design 
calculations that include the kinetic inductance of the 
Re center conductor. A spectroscopy measurement [3] 
yields a splitting size of SI /2/T = 19.0 ± 0.5 MHz, as 
expected for the qubit-resonator coupling strength due 
to a design value of coupling capacitor of 2.6 fF. The 
qubit-resonator coupling can be turned off by tuning the 
qubit off-resonance. In this experiment the detuning was 
400 MHz. 

The pulse sequence for generating and measuring Fock 
states [7j] is illustrated in Fig. [2^b). The timing of the 
pulses is based on the Jaynes-Cummings model [22j |. 
which predicts that on resonance, the occupation of the 
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FIG. 2: (a) Measurements of energy decay and Ramsey fringes 
for the resonator. Lines are fits to the data, yielding a res- 
onator T[ = 3.47 [is and T' 2 = 6.53 (b) Pulse sequence 
used for the experiment, showing both the flux bias (top) 
and microwaves (bottom) needed to generate the n = 5 Fock 
state, (c) Experimental data showing the state probability P e 
for the excited state |e) versus the qubit-resonator interaction 
time r during a Fock state readout (see text). The decrease 
in the qubit oscillation period from bottom to top, scaling as 
^Jn, is expected from the Jaynes-Cummings model for an n- 
photon Fock state. For each panel, the y-axis ranges from 0% 
to 85%. The purity of the sinusoidal oscillations indicates the 
high fidelity of the state preparation, (d) Fourier transform 
of the data in (c). For each Fock state, the Fourier amplitude 
(gray scale) is normalized to the maximum peak amplitude. 



two coupled degenerate states \g,n) and \e,n — 1), with 
n and n — 1 photons in the resonator, respectively, will 
oscillate with a radial frequency f2„ = i/nQ. For gen- 
eration, the n th photon is loaded from the excited state 
of the qubit using a carefully controlled interaction time 
of n/fl n — ir/^/nQ. The subsequent readout of the res- 
onator state is performed by tuning the qubit, initially in 
its ground state, into resonance with the resonator, and 
measuring the subsequent evolution of the qubit excited 
state probability P e as a function of the interaction time 
r. A Fourier transform (FT) of P c {t) gives the occupa- 
tion probability of the rt-photon Fock state as the Fourier 
component (in amplitude) of P e at the frequency Q n . 

Figure[2fc) shows P e versus interaction time for a num- 
ber of different Fock states. For these data, the readout 
starts after a delay time t ( i = 150 ns. For an interaction 
time r from to 300 ns, sinusoidal oscillations are clearly 
visible for Fock states from n=lupton = 15. The vis- 
ibility of the n = 1 sinusoidal oscillation is about 75%, 
which results from the high fidelity of the state prepa- 
ration and the high visibility of qubit readout close to 
90%. The frequency increases with increasing n, consis- 
tent with the predicted ^Jn dependence. A FT of the 
data is shown in Fig. [2£d), indicating the spectral purity 
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FIG. 3: (a) Qubit |e) state probability P e (gray scale) as a 
function of the qubit-resonator interaction time r (horizontal 
axis) and the measure delay td (vertical axis), (b) Fourier 
transform (along horizontal axis) of data in (a) , showing decay 
of the Fock states with increasing time. The Fourier spectrum 
is normalized to the maximum peak amplitude corresponding 
to n = 5 at td — 0. 
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of the oscillations. The Fock state n = 1 gives a Rabi 
frequency of 19.33 MHz, consistent with the splitting ob- 
tained from spectroscopy. 

The decay dynamics of the Fock states can be probed 
by varying the measurement delay td in the pulse se- 
quence. An example is shown in Fig. [3Ja), where P e is 
plotted in gray scale versus td and r for the Fock state 
n = 5. At td = 0, a periodic oscillation of P e is observed 
that corresponds to the occupation of the n — 5 Fock 
state. With increasing td, the oscillation evolves into a 
complex aperiodic pattern, finally changing to a slower 
oscillation frequency for td ~ 7 /is. The corresponding 
FT analysis is shown in Fig. [3jb). The system starts 
in the Fock state n = 5, and then evolves to a mixed 
state with several frequencies present in the spectrum, 
and eventually decays to a mixture of the Fock states 
n = 1 and n = 0. 

As no phase information is involved, the decay of Fock 
states can be simply described by changes in the state 
probabilities. The master equation is given by 

dP n (t) _ P n (t) | P n+1 (t) 
dt T n T n+ i 

where P n (t) and T n are the occupation probability and 
the decay time of the n-photon Fock state, respectively. 
To compare our data with theory, we plot in Fig. 0|e) 
the normalized amplitude in the FT as a function of the 
measure delay td for each Fock state in the spectrum of 
Fig. [3^b). Also shown in Fig. [4] is the decay of Fock 
states from n = 1 (a) to n = 6 (f), similarly obtained 
from the FT spectrum. The predictions from the mas- 
ter equation (lines) fit the data remarkably well for Fock 
states up to n = 6, taking as initial conditions the mea- 
sured state probabilities. The fitting parameters for T n 
for each panel are listed in Table |TJ The average val- 
ues show good scaling with the predicted T n ss T\ jn for 
n = 1 to 5. Since 600 ns-long time traces are used for 
the FT analysis to retrieve the Fock state probability, 
it is reasonable to begin to see some deviation from the 



FIG. 4: Time decay of Fock states from n = 1 (a) to n = 
6 (f), taken from the amplitude of the FT data. Occupation 
probabilities of the Fock states, P n , are plotted with symbols 
defined at top panel. Data are normalized by choosing the 
total occupation probability, including n — 0, at the measure 
delay td — to be unity [7|. Lines are fits to the data; the 
n-state lifetimes T„ as extracted from the fits are listed in 
Table U 

prediction for Tq w 640 ns. 

The FT analysis can also be used to test the decay dy- 
namics of a coherent state. A coherent state is a super- 
position of Fock states, with the occupation probability 
P n of the Fock state n given by a Poisson distribution, 

, \ a n e~ a 

P n (a) = j-, 2 

n ! 

where a is the average photon number. The resonator is 
first driven with a 100 ns classical microwave pulse, cre- 



TABLE I: Lifetimes (in units of fjs) of the n-photon Fock 
states taken from fits to Fig. UJa)-(f), respectively. The av- 
erages are listed in the second-to-last row, along with the 
expected T\/n scaled from the measurement of T\ = 3.20 fis. 
The values in parentheses are not included in the averages, 
as the fits are only reliable if there is a sizeable exponential 
decay of the data at times less than 7 /is, a limit set by our 
electronics. 
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0.80 
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a statistical mixture from a pure coherent state. A com- 
plete tomographic measurement is necessary in order to 
verify the phase coherence of the coherent state. 

In conclusion, we have improved the decay time of our 
earlier resonator design by a factor of 3, and thereby gen- 
erated high-fidelity Fock states with up to 15 photons. 
The time decay of both Fock and coherent states were di- 
rectly measured, and shown to be in excellent agreement 
with the theoretical prediction of a master equation with 
T n =Ti/n. 
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FIG. 5: (a) Occupation probabilities for coherent state analy- 
sis, taken from the normalized Fourier amplitude versus Fock 
state number n at a few example delay times. Data are nor- 
malized by choosing the total occupation probability at the 
measure delay td — to be unity, assuming that the n — 
probability at td — is negligible. Lines are fits to the data 
according to the Poisson distribution, (b) Average photon 
number a versus the measure delay td, obtained from fits to 
data as shown in (a). The line is an exponential fit yielding 
a lifetime of 3.43 /xs. 



ating a coherent state with a proportional to the drive 
pulse amplitude @. After a measure delay td, the res- 
onator state is read out using the qubit, initially in its 
ground state, in the same manner as discussed above. 
For a coherent state with (n) = a « 7.14, the occupation 
probabilities of the Fock states versus the photon number 
n are plotted in Fig.EJa) for several time delays. At each 
delay time td, the Fock state probabilities (points) can be 
well described by the Poisson distribution given by Eq. [2] 
(lines), with a as the fitting parameter. In Fig. G^b), we 
plot a as a function of td obtained for the complete data 
set. The average photon number a decays with a de- 
pendence that is fit by an exponential with a decay time 
3.43 /is, in good agreement with the measured resonator 
Tx. 

We note that although the occupation probabilities of 
the Fock states still follow a Poisson distribution during 
the entire decay process, our analysis cannot distinguish 
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